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Abstract 

We study the interaction of Gaussian one- and two-photon pulses with a single two-level 
atom based on a one-dimensional model of pulse propagation to and from the atom. The 
characteristic time scale of the atomic response is the dipole relaxation time l/T. We there- 
fore compare the effect of the non-linear two-photon interaction for a long pulse length of 
10/r with a short pulse of 1 /F. Our results indicate that the effect of the non-lineai^ interac- 
tion is particularly strong for the short pulse length of 1 /F. 

Keywords: Full quantum theory of light-atom interaction, Cavity QED, Non-linear photon- 
photon interaction, two-photon wavefunction. Single atom nonlinearity. One-dimensional 
atom. 

1 INTRODUCTION 

Nonlinearities are important in optical systems since they enable us to manipulate light by light. 
Especially, nonlinearities sensitive to single photons are necessary for quantum optical applica- 
tions such as optical quantum computation lH), entanglement generation 0, and quantum non- 
demolition measurements f3]. In particular, strong nonlinearities are indispensable for quantum 
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computation, since conditional operations at the single photon level are required to realize ele- 
mentary gate operations. However, conventional optical media require very high photon densities 
to obtain significant nonlinear effects f4l|. Single atom cavity quantumelectrodynamics may offer 
a possible solution 0161. In particular, a cavity can focus all the light in an input beam on a single 
atom, effectively creating a situation where a single atom is coupled to a one-dimensional light 
field |l7|. Such a one-dimensional atom is an ideal device for the non-linear manipulation of few 
photon pulses. It is therefore of great interest to investigate the response of a one-dimensional 
atom to one- and two-photon input pulses. 

There are many theories for such atom-cavity systems 0. However these theories usually 
eliminate the quantum state of the field outside the atom-cavity system and therefore the anal- 
ysis of the precise spatiotemporal coherence of the input and output photon pulses can not be 
achieved. In our previous work, we have therefore presented a theory for one- and two-photon 
pulses that includes the propagation to and from the system in a bad-cavity regime ||H1I2||, based 
on a one-dimensional model of light-atom interaction @. In these studies, we focused on the 
changes of the photon-photon correlations in time [Sj and in frequency 

In this paper, we study the effect of the non-linear interaction with the one-dimensional atom 
on the shape of one- and two-photon wavefunctions using Gaussian input pulses of different 
pulse lengths. For this purpose, we first introduce an analysis of the interaction in terms of 
absorption and transmission processes in section El The characteristic time scale defining the 
response of the atom is given by the dipole relaxation time 1 /V. In section|3l we then analyze the 
interaction processes for a long pulse length of lO/F in terms of pulse amplitude and pulse delay 
time. From these results, we can estimate that a maximal non-linear effect should be obtained 
for a pulse length equal to the dipole relaxation time of the atom. In section |4l we therefore 
analyze the corresponding interaction processes for a short pulse length of l/F, confirming the 
dominance of the non-linear effect throughout the short two-photon pulse. The results of the 
analysis are summarized and discussed in section|5] 

2 MODEL AND THEORY 

2.1 Theoretical model and physical realization 

The model of a single two-level atom in one-dimensional free space is shown in figure. [TJ a). In 
this model, an incoming one-photon or two-photon wavepacket interacts with the atom locally, 
and the output wavepacket is then emitted into the output field. The physical realization of this 
model can be implemented by a two-level atom coupled with a single mode of a one-sided cavity 
in the bad cavity regime QO. The cavity geometry is shown in figure. [T](b). The input field 
of the one-dimensional free space shown in figure, [l] (a) corresponds to the input of the one- 
sided cavity in figure. [T](b). Likewise, the output field in figure. [T](a) corresponds to the output 
of the one-sided cavity in figure. [11(b). Here, we assume that the spontaneous emission rate 
through the cavity mode is much larger than the spontaneous emission rate through non-cavity 
modes. In terms of the conventional cavity quantum electrodynamics parameters, this regime 
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is characterized by k ^ 5 ^ 7, where k is the cavity damping rate through the left mirror, g 
is the dipole coupling between the atom and the cavity mode and 7 is the rate of spontaneous 
emission into the non-cavity modes. Since the cavity damping rate k is much faster than the 
dipole coupling g, the method of adiabatic elimination can be applied to the time evolution of 
the cavity field ifTOll . This means that the interaction between the atom and the outside field 
through the cavity field can be expressed by an effective dipole relaxation rate T = g'^/n. The 
dipole relaxation rate T describes the dipole damping caused by emissions through the left mirror 
of the cavity, and the corresponding rate of spontaneous emission through the cavity is equal to 
2r ifTOlfTTl lTl. In our case, we assume that the rate of spontaneous emission into the non-cavity 
modes 7 is negligible (7 <^ 2r). Nearly all emissions from the atom can then be confined to 
the cavity and 2T is the total spontaneous emission rate of the excited atom in the cavity. In 
present cavity designs, this can be realized by exploiting the Purcell effect. For example, in the 
case of Turchette et al.'s experiment ["F], the cavity parameters indicate that about 80% of the 
sponaneous emission from the atom is emitted through the single cavity mode. 

(a) . (b) 

one- dimensional atom 



"a 
S 

<: 

^^^^^"^^ propagation 
input pulse 



two-level atom 




Figure 1: (a) Schematic of a two-level atom in one-dimensional free space and (b) physical 
realization of the one-dimensional atom using an atom-cavity system. An incoming one-photon 
or two-photon wavepacket interacts with the atom locally, and the output wavepacket is then 
emitted into the output field. The input field of the one-dimensional free space shown in (a) 
corresponds to the input of the one-sided cavity in (b). Likewise, the output field of the one- 
dimensional free space shown in (a) corresponds to the output of the one-sided cavity in (b). 



2.2 General solution of the interaction dynamics 

In the following, we investigate the effect of the one-dimensional atom on one and two-photon 
input pulses of different pulse length. For this purpose, we have to formulate the quantum states 
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for a given input pulse shape for both one photon and two photon cases. If the real space repre- 
sentation for the one-photon state reads, 



|A) = J Ma(x)|x), (1) 
the two-photon state can be expressed as the product state 

\ij)m = |Ai) ® IA2) = / c/xirfx2^'A(a;i) ■ ^'A(a;2)|a;i;x2). (2) 



The light pulses characterized by the above input states are propagating at the velocity of light 
in the input field of the schematic shown in figure. CJa). To simplify the expression, the spatial 
coordinates, Xi and X2, therefore represent a coordinate system moving at the velocity of light. 
The indices 1 and 2 of xi and X2 in eq. Q identify the two particles. The spatial features of the 
probability amplitudes given by the wavefunction correspond to the spatial features of the 
pulse mode A. Thus, the wavefunction ^a(2^i) • ^a(2^2) characterizes the single mode two-photon 
wavefunction, where the wavefunctions of both particles overlap perfectly. In our previous work 
lIHllZI, we derived the unitary time evolution operator in the Hilbert spaces for the one photon 
component and the two photon component by solving the Schrodinger equation of the interaction 
between a one-dimensional field and a single atom. The output one-photon state in the far field 
can then be obtained by convoluting the one-photon input wavefunction "^Aix) with the matrix 
description of the unitary operator in real space, 

|^)lphoton = j dx '^out{x)\x) 

with ^out(a;) = j dx' uiphoton(a;; a;') ■ ^a(x'). (3) 
The matrix elements uiphoton(a;, x') are given by 

2r r ' 

Uiphotonfa^; x') = S(x — x) e~~'^^ for x < x , else 0. (4) 

c 

Likewise, the output two-photon state in the far field can be obtained by a linear convolution, 

dXidX2^ oyA{.Xi] X2)\Xi] X2) , 



wi,h 4„.,fe;x.) = /dx;d.iu(x,x,;.;.4)-*.(x;).*.(4)^ (5) 
The matrix elements u(xi, X2; x'l, x'2) are given by 

u(xi, X2', Xi, X2) = Uiphoton(a^i; ^i) ■ Uiphoton(2;2; ^2) + Au^°"''"(xi, X2', Xi, X2), 
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where the contribution of the non-linear photon-photon interaction is given by 



Au^°"^'"(xi, X2] x'l, x^) = _i£!e-7(-'i +4-1-2) foj. x^^x2< Min[x;, x^]. (6) 



c2 



The output wavefunctions ^out(2;) and ^out(a;i; 2:2) describe the spatial feature of the output states 
in the far field of the atom shown in figure. [T] (a). Note that the spatiotemporal features of the 
input one-photon and two-photon wavefunctions are generally not preserved in the output due to 
the atomic response described by the matrix elements uiphoton(a;; x) and Au^°"^'"(a:i, X2] x'-^^x^). 

2.3 Components of the interaction 

The matrix element of the time evolution uiphoton(a;; x) in eq.© can be expanded in terms of 
two interaction processes, 

U-lphoton (-^j ) U-pi-op(x, X ) -|- Uabs('^) X ), (7) 

where the two processes are (I) single photon transmission without absorption, given by 

Uprop(a;; x) = 5{x - x), (8) 
and (II) single photon reemission after absorption, given by 

2r _r, 



e"^(^' ""^^ for x<x , else 0. (9) 



c 



The output one-photon wavefunction \l'out(2^) can thus be interpreted as the result of quantum 
interference of the two processes, 

^out(x) = ^prop(x) +^abs(x), (10) 

where the components are given by 

/oo 
rfx'up,op(x,x') ■ '^f^{x') = ^A(a;) 
-oo 

/oo 
rfx'uabs(x,x') ■ ^A(a;')- (11) 
-oo 

In the same way, the matrix element of the time evolution u(xi, X2; x'^, x'2) can be expanded 
in terms of three interaction processes. Process (1) is the transmission of both photons without 
absorption. Process (2) is the transmission of one photon without absorption and the absorption 
and reemission of the other photon. Process (3) is the absorption and reemission of both photons. 
In the following, we will refer to process (1) as the transmission process, to process (2) as the 
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one photon absorption process, and to process (3) as the two photon absorption process. The 
transformation of the two-photon wavefunction can then be written as 



u(a;i,a;2;a:^,x'2) = u'^^^(xi, X2; x^, X2) + u^^^(a;i, X2; x'^, X2) + u^^^(xi, 0:2; x'-^, X2). (12) 
The components of u(xi, X2; x^^x^) read 



2) 



U (Xi, X2, ^2) — Up].op(xi, X]^) ■ Upj.op(x2, X 

u^^^(a;i, a;2; x'^, x'2) = Uprop(xi; x[) ■ Uabs(x2; x'2) + Uabs(xi; x[) ■ Uprop(x2; X2) 
u^^^(xi, X2; x[, x'2) = Uabs(xi; x[) ■ Uabs(x2; X2) + Au^''"''" (xi , X2; x'^, X2). (13) 

Note that the non-linear contribution Au^°"''" only occurs in the two photon absorption process. 
This is because the non-linearity is due to the saturation of the two-level atom, which prevents 
the absorption of two photons at the same time. 

Likewise, the output wavefunction \l/out(xi, X2) can be expanded as 

^out(Xl,X2) = ^|,ut(Xi,X2) + ^|,'i(Xi,X2) + ^|,u't(a;i,X2), (14) 

where the \l/^'^(xi, X2) are given by 

/oo 
(ix'-LC/X2U^'^(Xi,X2,x'i,X2) ■ ^in(Xi,X2) (z=l,2,3). (15) 

-oo 

This analysis will be convenient to understand the origin of the statistical properties of the output 
photons. In the next section, these results are applied to the cases of Gaussian one-photon and 
two-photon input wavepackets. 



2.4 Definition of input pulses 

The dependence of the coherent spatiotemporal correlations between the output photons on the 
input pulse length is analyzed numerically using the results in the previous subsection. For our 
analysis, we use Gaussian one-photon and two-photon input wavefunctions, 

with = v/(7r/2)(cr)2. (16) 

The pulse length T is defined as twice the dispersion of the Gaussian distribution of the photons 
in time. Note that the response of the atom is governed by the time scale 1 /F, therefore the atom 
responds differently to long pulses and short pulses. In the following, the term 'long pulse' refers 
to a pulse whose pulse length T is much longer than the dipole relaxation time 1 /F. Likewise, the 
term 'short pulse' refers to a pulse whose pulse length T is about equal to the dipole relaxation 
time 1/F. In the next section, we analyze the output wavefunctions for long one-photon and 
two-photon pulses in order to identify the characteristic quantum interference effects. 
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3 RESPONSE TO A LONG INPUT PULSE 



3.1 Explanation of the one-photon output wavepacket in terms of one pho- 
ton processes 

Figure. 121(a) shows an example of a one-photon long pulse (broken line) and the corresponding 
output wavefunction (solid line). In this example, we have chosen an input pulse length of 
T = 10/r which is 10 times greater than the dipole relaxation time 1 /V. 

(b) „:-t- ' ., 




l-H 

o 
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Figure 2: Gaussian input and output one-photon wave packets for an input wave packet length 
of T = 10/F. The horizontal axes x represents the spatial coordinate of one photon. In (a), 
the broken line corresponds to the input wavefunction and the solid line corresponds to 

the output wavefunction ^out(2;)- In (b), the dotted line corresponds to the component of the 
process (I): single photon transmission without absorption and the thin line corresponds to the 
component of the process (II): single photon reemission after absorption. 

The pulse shape of the output one-photon pulse is almost identical with the input pulse ex- 
cept for the TT phase flip in the output amplitudes compared to the input amplitudes, '^oat{'^) ~ 
— ^A(a^)- Furthermore the minimum of the output amplitude is delayed compared to the max- 
imum of the input amplitudes. The delay time is about 2/F. Let us now analyze the origin of 
these features. 

Figure. |2l (b) shows the components corresponding to different light-atom interaction pro- 
cesses as given by eqs. dTTT) . The dotted line is the component of single photon transmission 
without absorption, ^! 



prop 



x) = ^!^{x) (see eqs. dTTT) ). The thin line is the component of single 
photon absorption, \E'abs(a^)- Since the absorption process convolutes the input amplitude with 
an exponential function according to eqs. (fTTT) . this component is not a precise Gaussian any- 
more. However, in the long pulse limit, the convolution function is so much narrower than the 
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Gaussian that the change in pulse shape is negligibly small. The component \l/abs(a;) can then be 
approximated as 

^abs(x) ^ -2^a{x) (17) 

As this approximation indicates, the amplitude of the absorption process is about two times as 
large as the amplitude of the transmission process. It is then possible to understand the vr phase 
flip in the output as a result of interference, 

^out(x) ^ ^a{xI -2^a(xH -^a(x), (18) 

transmission absorption 

where the negative amplitude of the absorption process is twice as high as the positive amplitude 
of the transmission process. 

However this approximation cannot explain the delay of the peak of the output wavefunction 
in figure. 2 (a). In the long pulse limit, the delay time of 1 /T in \I/abs(a^) can be included in the 
approximation given in eq.dTTt by expanding \l/abs(a^) as follows, 

^abs(x) ^ -2^A(x + c/r) 

~ -2^A(x)-2c/r^^A(x). (19) 
The delay time of \E'out(a^) can then be obtained from the interference of ^pmpix) and "^ahsix), 

*out(x) ^ *a(x) + (^-2^a{x) - 2c/r^$A(x; 

*a(x) + 2c/r^*A(x: 
^ -*A(a; + 2c/r). (20) 
The delay time of the output wavefunction is therefore equal to 2/T, as seen in figure. |2l(a). 



3.2 Explanation of the two-photon output wavepacket in terms of two pho- 
ton processes 

Let us now consider the output wavefunction for the long two photon pulse. The output wave- 
function can be analyzed in the same way as the one photon case, using the three interaction 
processes discussed in section |2l3. However, there are now two photon coordinates, xi and X2, 
and, as shall be seen in the following, the distance | Xi — X2 I between the two photons signifi- 
cantly influences the two photon absorption process. Note also that the bosonic symmetry of the 
indistinguishable photons requires that X2) = \&(x2, xi). 

Figure. |3] (a) shows the contour plot of a long two-photon input pulse \E'a(xi) ■ '^a{x2)- 
The probability amplitude increases from black to white shading. As before, the input pulse 
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has a length of T = lO/F. Figure. |3] (b) shows the cross-section of the input wavefunction at 
\xi — X2I = cr for time differences of r = 0, r = 1.4/r, andr = 5/r between the two photons. 
Note that the horizontal axis gives the average position of the two photons. The position of the 
peaks then defines the average delay of both photons at a given distance of cr between the two 
photons. 

Figure. |3](c) and figure. |3l(d) show the output wavefunction ^outl^^i, 2:2). A remarkable new 
feature of this output wavefunction is the negative amplitude region around xi = X2- The output 
amplitude changes from negative to positive values at about |a;i — X2\ = lAc/V. When going 
further away from X\ — X2, the spatial features of the output wavefunction eventually become 
similar to that of the input wavefunction. 

In order to illustrate these spatial features in detail, figure. |3l(d) shows the cross-sections of 
the output wavefunction at \xi — X2I = cr for r = r = 0, r = 1.4/r, and r = 5/r. The 
negative amplitude of the cross section at a time difference of r = (solid line) is about 3 times 
as large as the amplitude of the input wavefunction at r = shown in figure. |3](b) (solid line). 
Moreover, the peak of the output at r = is delayed by about 2/(3r) compared to the peak 
of the input wavefunction at the same time difference. The cross-section at a time difference of 
r = 1.4/r (broken line) is nearly zero. This means that the probability density of detecting one 
photon at a time difference of r = 1.4/r after detecting the other photon is vani shingly small at 
all times. The cross-section at r = 5/r (dotted line) is similar to the cross-section of the input 
wavefunction at the same time difference between the photons. However the output at r = 5/r 
is delayed compared to the input wavefunction by about 2/r, which is equal to the delay time 
observed in the one photon case. 

Let us now consider why such features appear. The two photon case can be understood as 
a superposition of three processes in the same way that the one photon case can be understood 
as a superposition of two processes. FigurelH shows the contour plots of the three components 
given by eq.(fT5l. Figure. El(a) corresponds to the transmission process, figure. 111(b) corresponds 
to the one photon absorption process, and figure. |4l(c) corresponds to the two photon absorption 
process, respectively. The component of the transmission process (figure. |4l(a)) is identical with 
the Gaussian input wavefunction, just as in the one photon case, ^out(a;i, X2) = ^a(2^i) ■ ^a(2^2)- 
Likewise, the component of the one photon absorption process (figure. |3| (b)), although not a 
precise Gaussian, can be approximated by the products of single photon transmission and single 
photon absorption. 



This approximation indicates that the negative amplitude of the one photon absorption process 
in figure. HI (b) is about 4 times as large as the positive amplitude of the transmission process 
in figure. 111(a). Moreover, the wavepacket of the one photon absorption process is delayed by 
c/ (2r) compared to the transmission process. 




{Xi,X2) 



(-2^a(xi + c/r)) • ^^{X2) + ^a(xi) ■ (-2^a(x2 + c/r)) 




A^a{xi + c/2r) ■ ^a(x2 + c/2r) 



(21) 
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Figure 3: Input and output wavefunction for a two-photon input pulse with a pulse length of 
T = 10/r. (a) is a contour plot of the input wavefunction \1'a(xi) ■ \1/a(x2). (b) is the cross 
section of the contour plot at \xi — X2I = cr, where the time differences are r = (solid line), 
T = 1.4/r (broken line) and r = 5/T (dotted line). The horizontal axis (xi + X2)/2 in (b) 
represents the average position of the two photons. Likewise, (c) is a contour plot of the output 
wavefunction ^^outl^^i, ^2) and (d) shows cross sections of the contour plot at \xi — X2I = cr, 
where the time differences are r = (solid line), r = 1.4/r (broken line), and r = 5/T (dotted 
line). 10 



The component of the two photon absorption process (figure. |4|(c)) includes the non-linear 
interaction between the two photons. For |xi — X2I ^ c/F, it resembles the shifted Gaussian 
expected for separate one photon absorption processes (see eq. (fT9t ). 

¥'^\xuX2) ^ 4^A(xi + c/F)-^A(a;2 + c/F). (22) 

However, the saturation of the atom completely suppresses this process for \xi — X2I = 0. Closer 
inspection shows that the transition between |xi — a;2| =0 and |xi — X2I ^ c/F can be approx- 
imated by an exponential function of |a;i — X2I, 

*ouU^i,^2) ~ 4^A(a:i + c/F)-*A(a:2 + c/F)(l-e-^l"^-"^l) (23) 

The combination of the two photon absorption process is therefore directly dependent on the 
time difference r = |xi — X2I /c between the two photons. 

In order to illustrate the dependence of the output component on the relative distance be- 
tween two photons in detail, figure. |5l shows the cross-sections at |xi — X2I = cr of the output 
amplitude \l/out(a;i, X2) (solid lines) and the components of the transmission process '$g^i{xi,X2) 
(dotted lines), the one photon absorption process '$''g^i{xi, X2) (thin lines), and the two photon 
absorption process ^'ouU^^ij 2:2) (broken lines). The time differences are r = in (a), r = 
1.4/F in (b), and r = 5/F in (c). For all three time differences r, the amplitude of the compo- 
nent of the one photon absorption process (thin lines) is about 4 times as large as the amplitude of 
the corresponding component of the transmission process (dotted lines). On the other hand, the 
ratio of the amplitude of the component of the two photon absorption process (broken lines) to 
the amplitude of the corresponding component of the transmission process (dotted lines) changes 
depending on the time difference r between the two photons as discussed in the following. 

Figure. |5] (a) shows the cross sections of the output components at r = 0. In this case, 
the amplitude of the component of the two photon absorption process, \E'out(xi, X2), is exactly 
zero (see eq. (l23t '). The output amplitude is then given by the sum of the transmission ampli- 
tude \I/out(xi, X2) = "^a{xi) ■ '^a{x2) and the one photon absorption amplitude \E'out(a;i, X2) ~ 
-4^'A(a;i) ■ ^^A(a;2), 

^out(a:i,X2) ^ $A(a;i) ■ *A(a;2) - 4^'A(a;i) ■ ^'a(x2) 

^ v ' ^ V ' 

transmission one photon absorption 

= -3^a(xi) ■ ^a(x2). (24) 

The amplitude of the output at the time difference of r = (solid line in figure. |5](a)) is therefore 
about 3 times as large as the corresponding amplitude of the input. 

Figure. |5](b) shows the cross sections of the output components at r = 1.4/r. At this special 

delay time, the factor of ^1 — e~ c l^i~^2l j describing the effect of the nonlinearity in eq. (l23t is 

approximately equal to 3/4. Therefore, the two photon absorption component in the total output 
can be approximated by '^^^\xi, X2) ~ 3'^a{xi) ■ '^a{x2). The output wavefunction is then given 
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Figure 4: Contour plots of the two photon processes W^Ji ° (xi, X2). (a) corresponds to the 
transmission process, (b) corresponds to the one photon absorption process, (c) corresponds to 
the two photon absorption process. 
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(a) (b) 




-20 -10 10 20 

{xi + X2)/2 (in units of c/r) 

Figure 5: Cross-sections of the output amplitude \1/ out (2^1, 3^2) and its components 
^out '° ^^\xi,X2) at time differences r = \xi — X2I /c of (a) r = 0, (b) r = lA/T and (c) 
r = 5/r. The solid lines correspond to the total output amplitude ^out(2;i, 2^2), the dotted lines 
correspond to the transmission amplitude ^out(a;i, X2), the thin lines correspond to the one photon 
absorption amplitude \I/out(a;i, X2), and the broken lines correspond to the two photon absorption 
amplitude \I/out(a^i, X2). 
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by 

^out(Xl,X2) ^ ^A(a:i) ■ ^a(2;2) -4^A(Xl)-^A(x2) + 3^A(a;i)-^A(x2) ^0. (25) 
^ V ' ^ V ' ^ V ' 

transmission one plioton absorption two photon absorption 

At a delay time of r = 1.4/r, the positive amplitudes of the transmission process and the two 
photon absorption process thus compensate the negative amplitude of the one photon absorption 
process, resulting in a total output amplitude of nearly zero. 

Figure. |5](c) shows the cross sections of the output components at r = 5/r. In this case, the 
amplitude of the component of the two photon absorption process, \l/out(xi, 0:2), can be approxi- 
mated as 4\E'A(a^i) ■ ^a(2^2) since the time difference r is sufficiently larger than 1 /F. The output 
amplitude is then given by 

^out(xi, X2) ^ ^a(xi) ■ ^a(x2) - 4^a(xi) ■ ^A(a;2) + 4^a(xi) ■ *A(a;2) 

V , ' -> ^ ' ^ ' 

transmission one photon absorption two photon absorption 

^ ^A(a;i) ■ *a(x2). (26) 

The amplitude of the output at the time difference of r = 5/F (solid line in figure. |5](c)) is thus 
almost identical with the amplitude of the input at the same time difference. 

To explain the delay of the output wavefunction, we can include the delay times in the argu- 
ments of the approximations given by eq. d^TT) and eq. (l23b . For r = 0, the output wavefunction 
can then be approximated as 

^out(a;i, X2) ^ ^a(xi) ■ *A(a;2) - 4^A(a;i + c/2r) ■ ^a(x2 + c/2r) 

^ v ' " V ' 

transmission one photon absorption 

/ 1c ( d d 

^ -3 ( ^a(xi) ■ ^A(a;2) + ^ ( ^^a(xi) ■ ^A(a;2) + ^a(xi) ■ — ^a(x2; 

^ -3^A(xi + 2c/(3r))-*A(x2 + 2c/(3r)). (27) 

The amplitude of the output wavefunction at a time difference of r = is thus delayed by 
about 2/(3r) compared to the amplitude of the input wavefunction at the same time difference. 
In the case of r = 1.4/r, the output wavefunction is nearly zero, and no delay time can be 
defined. For r = S/F, the two photon component ^'ouU^ii^a) of the output wavefunction can 
be approximated using the limit of \xi — X2I ^ c/F given in eq. (E^ . Since in this limit, the two 
photons are absorbed and reemitted independently of each other, the output wavefunction can 
then be factorized into one photon processes (see eq.lfTSt'). 

*out(a;i,X2) ^ (*a(xi) - 2*A(a;i + c/F)) ■ ($a(x2) - 2*a(x2 + c/F)) 

^ (-^a(xi + 2c/F)) ■ (-^A(a;2 + 2c/F)) . (28) 

The delay of the peak of the output wavefunction for \xi — X2I ^ c/F is thus equal to the 
delay for the one photon output wavefunction discussed in section |3ll. In summary, the delay 
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of the negative amplitude at xi = X2 is equal to 1/3 of the delay of the positive amplitude at 
\xi — X2I ^ c/r. This means that the effect of the non-linear interaction between the photons at 
xi = X2 reduces the total delay of the output wavefunction as compared to the delay caused by 
the uncorrected transmission and absorption processes at \xi — X2I ^ c/T. 



4 RESPONSE TO A SHORT PULSE 

As shown in the long pulse case, the atomic nonlinearity is characterized by a negative amplitude 
of the output wavefunction around xi = x^- In the long pulse case, this nonlinear effect is limited 
to a region of \x\ — X2I < 1.4c/r. In order to realize a non-linear phase change throughout the 
pulse, it is necessary to shorten the pulse to a pulse length comparable to 1 /F. In the following, 
we therefore discuss the response of the atom to a pulse of length T = 1 /F. 



4.1 Explanation of the output one-photon wavepacket in terms of one pho- 
ton processes 

The line plot in figure. |6l(a) shows the short one-photon pulse (broken line) and the corresponding 
output wavefunction (solid line). The length of the input wavepacket is T = 1/F. The pulse 

(a) i.5| I (b) l.5i I 




-4 -3-2-10 1 2 -4 -3-2-10 1 2 

X (in units of c/F) x (in units of c/F) 

Figure 6: Gaussian input and output one-photon wavepackets for an input wavepacket length of 
T = 1/F. In (a), the broken line corresponds to the input wavefunction pj^x) and the solid 
line corresponds to the output wavefunction \E'oMt(a^)- In (b), the dotted line corresponds to the 
component of the process (I): single photon transmission without absorption and the thin line 
corresponds to the component of the process (II): single photon reemission after absorption. 



shape of the output pulse is now quite different from the input pulse. Specifically, the output 
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amplitude is still positive for x > and only changes to negative values in the latter part of the 
pulse. Let us now analyze the origin of these features. 

Figure. |6l (b) shows the components corresponding to different light-atom interaction pro- 
cesses as given by eqs. lfTTT) . The dotted line is the component of single photon transmission 
without absorption, ^prop(2;)- The thin line is that of single photon reemission after absorption, 
^abs(a^)- As before, ^pmpix) = "^Aix) according to eqs. lfTTT) . On the other hand, the absorption 
component \[^abs(2^) (thin line) is given by the convolution with the exponential function of the 
dipole relaxation, given by eq.® and eqs. lfTTT) . The peak of this component is shifted by about 
0.6/r compared to the component of the transmission process. Furthermore, the peak of the 
amplitude is only about 1.3 times as large as the corresponding amplitude of the component of 
the transmission process. Therefore, the front part of the output wavefunction is dominated by 
the positive component of the transmission process \l/prop(a^), even though the major part of the 
output wavefunction is still dominated by the negative amplitude component of the absorption 
process \E'abs(3;). The large delay between the components and the reduced amplitude of the 
dominant negative absorption component thus cause the asymmetry of the output wavefunction 
seen in figure. |6l(a). 

4.2 Explanation of the two-photon output wavepacket in terms of two pho- 
ton processes 

Let us now consider the output wavefunction for the short two photon pulse. Figure. IT] (a) shows 
the contour plot of the short two-photon input pulse "^Aixi) ■ '^a{x2)- The probability amplitude 
increases from black to white shading. As before, the input pulse has a length of T = l/F. 
Figure. 111(b) shows the cross-section of the input wavefunction at \xi — X2I = cr for time 
differences of r = 0, r = O.S/F, and r = l/F between the two photons. Figure. IT] (c) and 
figure. |7](d) show the output wavefunction ^out(a;i, X2). The output amplitude is now negative 
everywhere. However, the pulse shape of the output shown in figure. |7](c) is quite different from 
the shape of the input pulse shown in figure. |7](a). In order to illustrate the spatial features of the 
output wavefunction in detail, figure. IT] (d) shows the cross-sections of the output wavefunction 
for time differences of r = 0, r = 0.3/F, and r = l/F between the two photons. Note 
that the shape of the cross-sections of the output in figure. |7] (d) is similar to that of the input 
cross-sections in figure. IT] (b). However, the amplitudes are quite different. Specifically, the 
output amplitude at a time difference of r = (figure. IT] (d), solid line) is larger than the input 
amplitude at the same time difference (figure. |7](b), solid line), while the output amplitude at a 
time difference of r = l/F (figure. IT] (d), dotted line) is smaller than the input amplitude at the 
same time difference (figure. IT] (b), dotted line). This indicates that the main difference in pulse 
shape between the input and output pulses is that the output amplitude decreases faster than the 
input amplitude between r = and r = c/F. On the other hand, the delay time of the peak of 
the output amplitude does not change much. Its value is about 0.4/F for all time differences r 
between the photons. 

Let us now consider why such features appear. Figure. [HI shows the contour plot of the three 
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Figure 7: Input and output wavefunction for a two-photon input pulse with a pulse length of 
T = 1/F. (a) is a contour plot of the input wavefunction \E'a(2;i) ■ ^a{x2)- (b) is the cross- 
section of the contour line plot at \xi — X2I = cr, where the time differences are r = (solid 
line), T = 0.3/F (broken line) and r = 1/F (dotted line). The horizontal axis {xi + X2)/2 
in (b) is the average position of the two photons. Likewise, (c) is a contour plot of the output 
wavefunction \E'out(2^i5 ^2) and (d) shows cross sections of the contour plot at |xi — X2I = cr, 
where the time differences are r = (solid line), r = 0.3/F (broken line) and r = 1/F (dotted 
line). 1'^ 



components given by eq. (fT5t . Figure. [8] (a) corresponds to the transmission process, figure. [8] 
(b) corresponds to the one photon absorption process and figure. [H (c) corresponds to the two 
photon absorption process. As before, the component of the transmission process (figure. [8] (a)) 
is identical with the input wavefunction, '^^^^^{xi, X2) = ^prop(a^i) ■ ^prop(2;2)- Likewise, the 
component of the one photon absorption process (figure. [8] (b)) can be described as 

^out(a;i,a;2) = ^prop(a;i) ■ «'abs(a;2) + ^^abs(a;i) ■ ^'prop(a;2) (29) 

Since the single photon transmission \l/prop(a^) is positive and the single photon absorption \[^abs(a^) 
is negative, the total amplitude of this component is negative for all xi, x^- 

The component of the two photon absorption process (figure. [81(c)) includes the non-linear 
interaction between the two photons. The positive amplitude of this component is therefore al- 
ways smaller than the product of the single photon absorption amplitudes, \l/abs(3^i) ■ ^abs(a^2)- 
Specifically, the two photon absorption component \l/out(xi, X2) is exactly zero at x\ = X2 and in- 
creases rapidly up to a maximal amplitude at about \xi — X2I = c/T. For larger photon distances, 
l^i — X2I > c/T, \l/out(a:i, 0:2) decreases as a consequence of the decrease in the amplitudes of the 
photon absorption processes, \E'abs(a^i) ■ ^abs(a^2)- Note however that, at its maximal amplitudes 
near \xi — X2I = c/T, the two photon absorption amplitude ^'ouU^^i, 2:2) is still significantly 
lower than the product of the single photon absorption amplitudes, \I/abs(2;i) ■ ^abs(2;2)- This 
means that the non-linear interaction of the photons significantly affects the whole output pulse 
in the case of an input pulse length of T = 1/T. 

In order to illustrate the dependence of the output component on the relative distance between 
two photons in detail, figure. |9lshows the cross-sections at |xi — X2 1 = cr of the output amplitude 
^out(a;i,X2) (solid line) and the components of the transmission process \l/out(xi, X2) (dotted 
line), the one photon absorption process '^''^^^{xi, X2) (thin line), and the two photon absorption 
process \l/out(xi, X2) (broken line). The time differences are r = in (a), r = 0.3/r in (b), and 
r = 1/r in (c). In all figures, the component of the one photon absorption process (thin line) is 
always negative and significantly greater than the component of the transmission process (dotted 
line). Furthermore, the ratio of the peak of the one photon absorption amplitude (thin line) to 
the peak of the corresponding transmission amplitude does not change very much depending on 
the time difference r between the two photons. On the other hand, the ratio of the peak of the 
two photon absorption amplitude (broken line) to the peak of the corresponding transmission 
amplitude (dotted line) changes significantly depending on the time difference r between the 
two photons as discussed in the following. 

Figure. |9l (a) shows the cross-section of the output components at r = 0. In this case, the 
two photon absorption amplitude (broken line) is exactly zero. The output amplitude is then 
given by the sum of the transmission amplitude (dotted line) and the one photon absorption 
amplitude (thin line). Since the one photon absorption amplitude is about 2.5 times as large as 
the transmission amplitude, the total output pulse (solid line) is negative and about 1.5 times as 
large as the transmission amplitude. 

Figure. |9l (b) shows the cross-sections of the output components at r = 0.3/r. The one 
photon absorption amplitude (thin line) and the transmission amplitude (dotted line) are not much 
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Figure 8: Contour plots of the two photon processes ^'ouV° '^^''{^ii X2)- (a) corresponds to the 
transmission process, (b) corresponds to the one photon absorption process, (c) corresponds to 
the two photon absorption process. 
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different from r = (figure. |9l(a)). However, the two photon absorption amplitude (broken line) 
is now about half as high as the transmission amplitude (dotted line). Therefore, the output 
amplitude is reduced to about the same height as the transmission amplitude. 

Figure. |9l(c) shows the cross-sections of the output components at r = 1 /V. The two photon 
absorption amplitude (broken line) and the transmission amplitude (dotted line) are now sig- 
nificantly lower than for r = (figure. |9l (a)). However, the negative one photon absorption 
amplitude is still about 3 times as large as the transmission amplitude. The two photon absorp- 
tion amplitude (broken line) has increased to about 1.5 times the transmission amplitude, leaving 
a total negative output amplitude (solid line) of about 0.5 times the transmission amplitude. Thus 
the output amplitude is reduced to very low values, even though the components themselves are 
still very large. The fast decrease of the total output amplitude shown in figure. |7] (c) and |7] 
(d) can then be understood as a consequence of the rapid increase of the two photon absorption 
amplitude with r. 

Interestingly, the delay times of the pulse components are not sufficiently high to cause the 
same change in the sign of the amplitude observed for the one photon case (figure. |^(a)). No part 
of the two-photon output wavefunction is dominated by either the transmission or the two photon 
absorption amplitude. This is mostly because the two photon absorption amplitude is suppressed 
by the non-linearity. In fact, it can be seen in figure. |9l(c) that the nonlinearity even reduces the 
delay of the two photon absorption amplitude from a delay time of 0.6/r for \E'abs(a^i) ■ ^abs(a;2) 
to only about 0.45/r for ^out(xi,X2) at r = l/F. Therefore, the change of the pulse shape in 
the two photon pulses is actually less drastic than the change of pulse shape observed in the one 
photon case. 

5 SUMMARY AND CONCLUSION 

We have studied the effect of the strong nonlinearity of a one-dimensional atom on one- and 
two-photon input wavepakets. The response of the atom can be characterized by transmission 
and absorption processes. For long pulses, the absorption components dominate in the out- 
put. However, in the two photon case, the two photon absorption component is suppressed 
around xi = X2 by the saturation of the two-level atom. The long two photon pulse there- 
fore has a non-linear region of negative amplitude dominated by the one photon absorption for 
\xi — X2I < 1.4c/r, and a linear region of positive amplitude dominated by the two photon 
absorption for \xi — X2I > lAc/T. For a short pulse length of T = l/F, the two photon ab- 
sorption is strongly suppressed throughout the whole output pulse. The one photon absorption 
process then dominates everywhere, and the amplitude of the short two-photon output pulse is 
completely negative. It is therefore possible to obtain a very strong non-linear effect by adjusting 
the pulse length T to the dipole relaxation time of the atom. 

Besides the total amplitude, we have also characterized the change in delay time and pulse 
shape. For long pulses, the delay of the non-linear region is reduced to 1/3 of the delay of 
the liner region. In the long pulse case, the interaction of the two photons thus significantly 
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Figure 9: Cross-sections of the output amplitude \1/ out (2^1, 2^2) and its components 
^out '° ''^■'(^153^2) at time differences r = \xi — X2I /c of (a) r = 0, (b) r = 0.3/r and (c) 
r = 1 /r. The solid lines correspond to the total output amplitude ^out(2;i, 2^2), the dotted lines 
correspond to the transmission amplitude ^out(a;i, X2), the thin lines correspond to the one photon 
absorption amplitude \I/out(a;i, X2), and the broken lines correspond to the two photon absorption 
amplitude \&out(3;i; 2:2) • 
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reduces the average delay time caused by the absorption and reemission processes at the atom 
and causes a faster reemission into the output field. In the case of a short one photon pulse, the 
delay between the transmission and the absorption component is so great that the pulse shape 
is changed significantly. In particular, the single photon output amplitude has both negative 
and positive regions. On the other hand, the output pulse shape for a short two photon pulse is 
changed much less due to the reduced delay caused by the non-linearity. It is therefore possible 
to obtain a negative output amplitude throughout the short output pulse. This result indicates that 
the atomic nonlinearity is particularly strong at a pulse length of T = 1 /V. 

The analysis of the output in terms of components corresponding to different interaction 
processes thus allow us to understand the origin of the non-linear photon-photon interaction in 
great detail. These insights may be very useful for the development of novel non-linear optical 
devices operating at the few photon level. 
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